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Energies, transition rates, and electron electric dipole moment (EDM) enhancement factors are 
calculated for low-lying states of Ce IV and Pr V using relativistic many-body perturbation theory. 
This study is related to recent investigations of the more complicated Gd IV ion, which is promising 
for electron EDM experiments. The ions Ce IV and Pr V both have a single valence electron, 
permitting one to carry out reliable ab-initio calculations of energy levels, transition rates and other 
atomic properties using well developed computational methods. 
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I. INTRODUCTION 



Recently, searches for an electron electric dipole mo- 
ment (EDM) in atoms and molecules have gained con- 
siderable interest. Since an electron EDM violates time- 
reversal symmetry and parity, the discovery of an EDM 
would have many implications for modern fundamen- 
tal theories (a general overview is given by Khriplovich 
and Lamoreaux Q). The best limit on the value of the 
electron EDM was obtained by Regan et al. Q in Tl, 
d e < 1.6 x 1CP 27 e cm. Although there is a large en- 
hancement (~ 600) of the electron EDM in Tl, the den- 
sity of atoms in a beam is much lower than in a solid- 
state system. To utilize the advantage of high densities, 
it has been proposed to use Gadolinium Gallium Gar- 
net Gd3GasOi2 or Gadolinium Iron Garnet GdaFesO^ 
(densities ~ 10 22 /cc) in EDM experiments || §]. There 
are various features of these materials that make them 
particularly useful for such experiments. For example, 
magneto-electric effects are forbidden owing to the FCC 
symmetry of the crystals, simplifying the exclusion of 
systematics; the crystal GdsFesOi^ has a very high re- 
sistivity (> io 16 n -cm); and spin alignment is relatively 
easy. The electron EDM enhancement for the Gd 3+ ion 
(^2 — 3) is much smaller than for Tl; nevertheless there 
is a substantial gain in the number of atoms. 

Motivated by proposed EDM measurements, calcula- 
tions of EDM enhancement have been performed recently 



id . e du / ~ i s a vuk 



TRlectror 



*F,lcctromir. addres s- isavukov@princeton.edu; URL: http://www 



johnson 



ic address: 



johnson@nd.edu; URL: littp : //www.nd. edu/ 



* Electronic address: jsafrono@nd.edu 



^Electronic address: nsafrono@nd.edu 



|| q|. Core polarization effects in Gd IV were discussed 
by Dzuba et al. @ and found to decrease the size of 
EDM enhancement factor from -3.3 to -2.2 . Thus, there 
is a strong sensitivity of the EDM enhancement to many- 
body effects and further investigation is necessary. 

The ion Gd +3 , which has a 4/ 7 ground-state configu- 
ration is extremely complicated and difficult for atomic 
theory, though some understanding was gained by Dzuba 
et al. @] using a relativistic configuration-interaction 
(RCI) method and the widely used Cowan || code. 
As mentioned above, core polarization was found to be 
substantial. For example, to match the experimental 
4/ — > 5d transition energies, scaling factors of 0.8 and 
0.85 for Coulomb integrals were used in the RCI and 
Cowan codes, respectively. In addition, a polarization 
potential was introduced in RCI code to match experi- 
mental energies and a systematic energy shift of ~ 18, 000 
cm -1 was made in the Cowan code. This shift can be at- 
tributed to core polarization by a single valence electron. 
An ab-initio investigation of correlation should help to 
understand these empirical adjustments. The simplest 
ions (those having a single 4/ valence electron) which 
exhibit core-polarization effects similar to Gd IV are Ce 
IV and Pr V. To gain a clearer understanding of the role 
of core polarization and other many-body effects in Ce IV 
and Pr V, we calculate energies, dipole transition matrix 
elements, and EDM enhancement factors for low lying 
states using relativistic many-body perturbation theory 
(MBPT). 

The present calculations of energies are carried out to 
third order in MBPT using the methods developed ear- 
lier to study the Li, Na, and Cu isoelectronic sequences 
H [U^ 0. The calculations of transition matrix ele- 
ments is also carried out to third order in MBPT using 
the methods developed in to treat the alkali atoms 
and alkalilike ions. Here, we use the gauge-independent 
version of the third-order MBPT code described in fn3|. 
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The calculations of electron EDM enhancement factors, 
which involve a sum over intermediate states, are carried 
out in the random-phase approximation (RPA) following 
a procedure similar to that described in [[yj . 

It should be mentioned that Ce +3 , which has a 4/5/2 
ground state and an observed magnetic moment 2.3-2.5 
fj,B can be embedded into garnet crystals and used in 
EDM experiments of the type proposed in B. In those 
experiments, a strong electric field polarizes ions which 
in turn produce a small magnetic field that is measured 
in a sensitive SQUID detector. The magnetic field at low 
temperatures is proportional to the product of the elec- 
tric dipole moment of the ion and its magnetic moment. 
The estimated EDM of Gd +3 as is a factor of about 3 
larger than for Ce +3 and the magnetic moment of Gd +3 
(7.9-8.0 hb) is also about three times larger; therefore, 
there is an overall advantage of nine in Gd +3 compared 
to Ce +3 . Nevertheless, the Ce +3 ion could still compete 
in setting experimental limits on the electron EDM owing 
to the fact that its ionic EDM has significantly smaller 
theoretical uncertainty. 



II. CALCULATION OF ENERGIES 

First-, second- and third-order Coulomb energies 
n — 1..3, and first- and second-order Breit ener- 
gies fiW, n = 1,2, calculated using methods described 
in Refs. ||, O, are presented in Table | along with 
the resultant theoretical energies E to t and predicted ener- 
gies -Enist from the National Institute of Standards and 
Technology (NIST) given by Martin et al. |]|. We see 
that second-order corrections are large and improve the 
accuracy of the first-order Dirac-Hartree-Fock (DHF) en- 
ergies. Third-order MBPT further improves the ground 
state energy. However, the third-order correction is rela- 
tively large (roughly one-third of the second order) and 
overshoots the experimental value, which indicates that 
oscillations in higher orders are likely. Using a geomet- 
ric progression with q = —1/3 we can extrapolate the 
second- and third-order values to give a limiting ground- 
state energy —296420 cm -1 for Cc +3 , in close agreement 
with the predicted value from flB]] . The corresponding ex- 
trapolation for Pr +4 gives —468290 cm -1 differing from 
the predicted ground-state energy given in fl5f| by about 
4000 cm . This rather large difference casts doubt on 
the threshold energy of Pr V predicted in jl5| and ex- 
plains the large differences with the NIST energies seen 
111 lower half of Table § 

In Refs. §, [k], 0, where the ionic ground-states were 
2s, 3s, and 4s, respectively, the third-order correction 
was uniformly much smaller than the second-order cor- 
rection. The relatively large size of the third-order cor- 
rections in Ce IV and Pr V arise because of the double- 
well feature of the 4/ Coulomb potential discussed, for 
example, by Cheng and Froese-Fischer . 



TABLE I: First-order (DHF) energies E m , second- and third- 
order Coulomb energies E^ and E 1 - 3 ^, first- and second-order 
Breit corrections B' 1 ^ and S' 2 ' and totals -Etot for Ce IV and 
Pr V are compared with predicted energies -Bnist given by 
Martin et al. Ol. Units: cm" 1 



nlj 


E m 




E (2) 


B (2) 


E (3) 


Etot 


£nist 


5E 


Ce IV 


4/5/2 


-261361 


616 


-46747 


-2552 


15672 


-294372 


-296470 


2098 


4/7/2 


-259378 


423 


-46091 


-2511 


15398 


-292158 


-294217 


2059 


5d 3/2 


-236793 


292 


-13516 


-545 


3753 


-246809 


-246733 


-76 


5d 5/2 


-234637 


220 


-12957 


-528 


3675 


-244226 


-244244 


18 


6S1/2 


-203245 


168 


-10201 


-200 


3680 


-209794 


-209868 


70 


6P1/2 


-168978 


167 


-7283 


-127 


2329 


-173891 


-173885 


-6 


6^3/2 


-164703 


120 


-6648 


-122 


2118 


-169236 


-169178 


-58 


6d 3 / 2 


-116268 


67 


-3989 


-106 


997 


-119300 


-119272 


-28 


6d 5 / 2 


-115581 


51 


-3922 


-106 


977 


-118580 


-117557 


-1024 


7Si/ 2 


-110878 


65 


-3601 


-74 


1293 


-113196 


-112968 


-226 



Pr V 



4/5/2 


-431686 


798 


-48799 


-3008 


16527 


-466167 


-464000 


-2167 


4/7/2 


-428863 


551 


-48175 


-2965 


16270 


-463182 


-460973 


-2209 


5d 3 / 2 


-341470 


395 


-14968 


-669 


4088 


-352624 


-348948 


-3676 


5^5/2 


-338314 


298 


-14386 


-651 


3827 


-349225 


-345486 


-3739 


6S1/2 


-281133 


225 


-11649 


-247 


4241 


-288563 


-285029 


-3535 


6fl/2 


-239154 


239 


-9136 


-169 


2451 


-245769 


-240522 


-5247 


6P3/2 


-233035 


171 


-8489 


-163 


2040 


-239477 


-233961 


-5516 


6d 3 / 2 


-171610 


98 


-5746 


-141 


792 


-176576 






6d 5 /2 


-170572 


75 


-5647 


-141 


771 


-175481 






7Si/ 2 


-156266 


92 


-4489 


-97 


1329 


-161431 


-159489 


-1942 



III. CALCULATION OF TRANSITION MATRIX 
ELEMENTS AND TRANSITION RATES 

Transition matrix elements provide another test of 
quality of atomic-structure calculations and another mea- 
sure of the size of correlation corrections. Third-order 
MBPT reduced matrix elements for transitions between 
low-lying states of Ce +3 and Pr +4 are presented in Ta- 
ble [n| The first-order reduced matrix elements Z^ are 
obtained from length-form DHF calculations. Length- 
form and velocity-form matrix elements differ typically 
by 10%. Second-order matrix elements in the table Z^ 2 \ 
which include Z^ x \ are extended to include all higher- 
order corrections associated with the random-phase ap- 
proximation. These second-order calculations are prac- 
tically gauge independent. In the present calculations, 
length- and velocity-form matrix elements in the RPA 
agree to six or more digits. The third-order matrix el- 
ements Z^ include Z^ plus Brueckner-orbital (BO), 
structural radiation, and normalization corrections de- 
scribed, for example, in [|l2|. These calculations are car- 
ried out in a gauge-independent manner, including ap- 
propriate derivative terms, as described in [ jDl , We trun- 
cated our basis set to include only those partial waves 
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TABLE II: Reduced matrix elements of the dipole operator in 
first-, second-, and third-order perturbation theory for tran- 
sitions in Ce IV and Pr V. 











Ce IV 






Pr V 




Transition 


Z (D 


Z (2) 


Z (3) 


Z (i) 


Z (2) 


Z (3) 


4/5/2 




5rf 3 /2 


1.498 


0.972 


1.172 


1.146 


0.678 


0.706 


4/5/2 




5rf 5 /2 


0.396 


0.264 


0.308 


0.302 


0.186 


0.189 


4/7/2 




5^5/2 


1.799 


1.193 


1.413 


1.370 


0.828 


0.852 


5d 3 / 2 




6P1/2 


1.976 


1.768 


1.682 


1.648 


1.471 


1.396 


5d 3 / 2 




6P3/2 


8.363 


0.756 


0.719 


0.690 


0.625 


0.592 


5^5/2 




6P3/2 


2.585 


2.352 


2.229 


2.140 


1.950 


1.837 


6S1/2 




6P1/2 


2.847 


2.482 


2.402 


2.560 


2.207 


2.136 


6S1/2 


— ► 


6P3/2 


4.012 


3.512 


3.401 


3.609 


3.125 


3.020 



with I < 8, and found that length- and velocity-form 
third-order reduced matrix elements agreed to 4 digits. 

As can be seen in Table |J, RPA corrections are very 
large, 10-40%, being largest for 4/ — > 5d transitions, and 
must be taken into account. Such behavior can be at- 
tributed to core shielding which is substantial because 
valence electrons penetrate deeply into the core. Third- 
order corrections are smaller, 2-4% scaling as l/Zi on . If 
such scaling holds in higher orders, we can estimate the 
accuracy of our calculations to be 0.4-0.8%. The domi- 
nant contribution in third order comes from the BO cor- 
rection which is approximately equal to the sum of the 
other third-order corrections. 

Transition rates A (s _1 ), oscillator strengths /, and 
wavelengths A (A) for electric dipole transitions between 
low-lying states of Ce IV and Pr V are given in Table H. 
These data are calculated using the dipole matrix ele- 
ments Z^ from Table || and predicted NIST transition 
energies fla] . In the two final columns of Table [II , we 



compare our MBPT wavelengths with the wavelengths 
from Ref. pUf . We also compare our MBPT oscillator 
strengths with theoretical oscillator strengths obtained 
by Migdalek and Wyrozumska p7j . The data in JlTj 
were obtained using a relativistic model potential (RMP) 
approach together with a core-polarizatio n ( CP) model 
potential. Our data and that from Ref. [[17) agree well 
for bd — 6p and 6s — 6p transitions but differ for 4/ — 5d 
transition where / values are very small. 



TABLE III: MBPT transition rates A (s _1 ), oscillator 
strengths /, and wavelengths A (A) for transitions in Ce IV 
and Pr V. MBPT (a) oscillator strengths are compared with 
theoretical calculations (b) performed in Ref. [Q. MBPT (a) 
wavelengths are compared with wavelengths (c) predicted by 
NIST pi. 



Transition 






f m 


A («0 


A (c) 


Ce IV 


4/5/2 




5d 3 / 2 


8. 56 [7] 


0.0346 


0.0159 


2102 


2011 


4/5/2 




5d 5 / 2 


4.56[6] 


0.0025 


0.0013 


1994 


1915 


4/7/2 




5d 5 /2 


8.41 [7] 


0.0379 


0.0186 


2086 


2001 


5d 3 / 2 




6P1/2 


1.1119] 


0.157 


0.158 


1371 


1373 


5^3/2 


— > 


6P3/2 


1.22 [8] 


0.0304 


0.0206 


1289 


1289 


5^5/2 


— > 


6^3/2 


1.06[9] 


0.189 


0.189 


1334 


1332 


6S1/2 




6P1/2 


2. 72 [8] 


0.315 


0.327 


2785 


2779 


6S1/2 


— > 


6p 3 /2 


3.95[8] 


0.715 


0.745 


2465 


2458 


Pr V 


4/5/2 




5^3/2 


3.85[8] 


0.0290 


0.0318 


881 


869 


4/5/2 




5d 5 / 2 


2.01[7] 


0.0021 


0.0019 


863 


844 


4/7/2 




5^5/2 


3. 78 [8] 


0.0318 


0.0285 


878 


866 


5^3/2 


— > 


6P1/2 


2. 52 [9] 


0.161 


0.165 


936 


922 


5d 3 / 2 




6P3/2 


2. 70 [8] 


0.0306 


0.0311 


884 


870 


5d 5 /2 


— > 


6P3/2 


2.37[9] 


0.191 


0.193 


911 


897 


6S1/2 


— > 


6P1/2 


4.07[8] 


0.308 


0.321 


2337 


2247 


6S1/2 




6P3/2 


6.15[9] 


0.707 


0.746 


2037 


1958 



many-electric atom D may be written [|14| 
(vj v \eZ\n) (n\H cd 



E v — E n 

where eZ is the dipole electric operator 
eZ = ^2 ez n 



(1) 



and i?cdm is an equivalent EDM interaction 1 14 given by 



dm * 



IV. EDM ENHANCEMENT 

A. Basic equations 

According to Schiff's theorem the electric dipole 
moment of an atom induced by an intrinsic electron EDM 
vanishes in the nonrelativistic limit; however, as shown 
by Sandars ||l9| , the atomic EDM is nonvanishing rela- 
tivistically and can be a large multiple of the intrinsic 
electron moment for heavy atoms. If we assume that the 
electron has an intrinsic EDM d e , then the EDM of a 



This equivalent interaction, which automatically ac- 
counts for Schiff's theorem, is rotationally invariant and 
therefore conserves angular momentum; it violates both 
parity and time-reversal symmetry. 

For an atom or ion with one valence electron, one- 
electron matrix elements of i? c dm may be written in low- 
est order as 

{n\\H c dm\\v), (2) 
where the (somewhat unconventional) reduced matrix el- 
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ement is 

(n||-ffedm||v) : 

2c- / dr 

e Jo 



G n (r) 



d 2 F v _ k v (k v - 1) 
dr 2 r 2 



G n (r)F v (r) 



+F n {r) 



d 2 G v k v (k v + 1) 
dr 2 r 2 



F n {r)G v {r) 



In the above equation, Gfc(r) and Fk(r) are the large 
and small components, respectively, of radial Dirac wave 
functions. Similarly, we may write 



(vj v \eZ\nj v ) 



Jv 



(2j v + l)(j„ + 1) 



(v\\eZ\\n) (3) 



with 



(v\\eZ\\n)=e(K v \\C 1 \\-K v )x 

O 

rdr[G v (r)G n (r)+F v (r)F n (r)], 



C\ q {r) being a normalized spherical harmonic. The ex- 
pression for the atomic dipole moment in lowest-order 
MBPT then reduces to 



D« = 2 



Jv \ - 

(2j„ + l)(j„ + l) ^ 



(v\\eZ\\i) (i\\H edm \\v) 



.(4) 

where e k are eigenvalues of the valence-electron Dirac 
equation. 



B. RPA correlation corrections 

Lowest-order calculations of the induced atomic EDM 
are carried out in a frozen-core T^ -1 DHF potential. 



14[ | to be very sensitive 
lat reason, the lowest- 



Such calculations were shown in 
to correlation corrections. For t : 
order "bare" matrix elements in Eq. (Q) are replaced by 
"dressed" RPA matrix elements. 



1. Z-RPA 

Thus, we replace the lowest-order dipole matrix ele- 
ment in Eq. (||) by 

(w\\Z™ A \\v) = {w\\Z\\v) 

1)n -n+il {a||^ RPA ||n)^i(^q) 

an 

, 1)8 -r 1+ ll M^g^ i (5) 

where the index a extends over all core orbitals and 
the index n extends over all virtual orbitals permitted 



by angular-momentum selection rules. The quantities 
Zj(ijkl) are Coulomb integrals 



Zj(ijki) = Xj(ijki) + [J]£{; J J L j x L (ijik) , 



where Xj(ijkl) are defined by 

Xj(ijkl) = (-l) J (i\\Cj\\k) (j\\Cj\\l) Rj{ijkl) 
Rj(ijkl) being a Slater integral 



(6) 



Rj(ijkl) 



J '< 



JO 



drdr 



,J+ 



T [ Gl (r)G k (r)+F l (r)F k (r)] 



[G j (r>)G l (t') + F j (r>)F l (T>)]. 



We designate the corresponding approximation to 
the atomic EDM by Z? RPA . Note that if we replace 
(n||Z RPA ||a) by (n||Z||a) on the right hand side of 
Eq. (||), then we obtain the second-order correlation cor- 
rection to the valence-excited dipole matrix element. The 
atomic EDM calculated in this approximation is desig- 
nated by D z . 



2. H-RPA 

Similarly, we replace the bare matrix element of the 
EDM interaction (w||i? c dm||w) by its dressed counterpart 

(w\\H^ A \\v^ {w\\H cAm \\v) 

Ua] (a\\H^\\n)Z Q (wnva) 




[jv] £a ~ £r, 

\j a ) Z (wavn)(n\\H^ A \\a) 



(7) 



We designate the approximation to D obtained using the 
dressed matrix element from Eq. ([?]) by Dff A . Again, 
if we replace dressed matrix elements by bare matrix el- 
ements on the right hand side of Eq. (0), we obtain a 
second-order approximation to (w\ |-ff c dm| \v}- The result- 

(2) 

ing correction to the atomic EDM is designated by D H . 
Core-excited matrix elements (n||Z RPA ||a) and 
) and (0) satisfy sets of coupled 
y m EM- 



rRPA 1 1 „\ in Eqs | 
equations given explicit 



(nllfl^lio) 



C. Calculations of EDM enhancement factors 

The sums over intermediate states in Eqs. (f|j7]) are 
carried out using basis functions obtained as linear com- 
binations of B-splines as described in We use 40 
splines of order 7 and constrain the ions to lie in a cavity 
of radius R = 35 a.u. for Ce +3 and 30 a.u. for Pr +4 . 

A detailed breakdown of the contributions to D for 4/ 



states of Ce IV is given in Table [V, where we list the 
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TABLE IV: Comparison of first-order, second-order, and RPA 
calculations of the atomic EDM enhancement factor D/d e for 
4/ states of Ce IV. 

state £>W AD™ ADf D™ A£)| PA APf PA D RPA 
4/5/2 -0.382 -0.388 0.332 -0.438 -0.785 0.387 -0.780 
4/7/2 -0.002 -0.033 0.022 -0.013 -0.045 0.015 -0.032 



TABLE V: EDM enhancement factors D/d e for low-lying 
states of Ce IV and Pr V. 



State 


£,(1) 


AD^j PA 


adF a 


d rpa 


Ce IV 


4/5/2 
4/7/2 
5d 3 / 2 
5d 5 / 2 

6S1/2 

6P1/2 

6P3/2 


-0.382 
-0.00225 
-1.95 
0.425 

120. 
-158. 

2.89 


-0.785 
-0.0451 
-3.38 
-0.347 
27.5 
-30.9 
7.24 


0.387 
0.0151 
0.779 
-0.136 
-19.8 
19.8 
-1.07 


-0.780 
-0.0323 

-4.55 
-0.0628 

128. 

-169. 

9.06 


Pr V 


4/5/2 
4/7/2 
5d 3 / 2 
5rf 5 / 2 

6S1/2 
6j3l/2 
6P3/2 


-0.142 
-0.00266 
-1.80 
0.174 

127. 
-157. 

2.78 


-0.0806 
-0.00444 
-2.91 
-0.827 
27.3 
-30.0 
6.06 


0.0926 
-0.000912 
0.655 
0.125 
-22.2 
21.7 
-1.01 


-0.130 
-0.00802 
-4.05 
-0.528 

132. 
-166. 

7.83 



DHF approximation, , the second-order correction 
D^ 2 \ the RPA approximation, D RPA , and the individual 



contributions to the second-order and RPA corrections 
from the dipole and weak-interaction matrix elements. 
One can see from the table that the correlation correc- 
tions to the weak-interaction matrix element are compa- 
rable to or larger than the lowest order matrix element. 
Moreover, there are significant changes in these correla- 
tion corrections going from second-order MBPT to full 
RPA calculations. 

Finally, in Table |V|, we present DHF and RPA values 
of the EDM enhancement factors D jd e for the low-lying 
4/, 5d, 6s, and 6p states of Ce IV and Pr V. 



V. CONCLUSION 

We have studied energies, transition probabilities, and 
EDM enhancement factors for Ce +3 and Pr +4 . We found 
that perturbation theory converges quite slowly and that 
RPA corrections are the dominant correlation correc- 
tions for transitions. We use our third-order MBPT with 
"dressed" matrix elements to obtain accurate transition 
rates. The most interesting discovery is that RPA cor- 
rections modify lowest-order values of the EDM enhance- 
ment factor significantly. 
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